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ABSTRACT 


folictememtereat model 15 developed to determine the heat 
Mmeamster rate Of a cylindrical condenser section of a rotating 
heat pipe. This model is coupled to an existing code and an 
analysis is accomplished on both a smooth and axially finned 
condenser. The results of this analysis are compared to those 
of a similar analysis performed on a tapered condenser heat 
pipe using identical geometric and operating parameters. 

Results of the comparison indicate cylindrical condensers 
miwenwa Lower heat transfer rate than equivalent tapered con- 
densers. This reduction in heat transfer rate is due to a 
greater condensate film thickness and is most significant in 
m@eomooth condenser. 

Axially finned condensers with triangular and rectangular 
fin profiles are also compared. The rectangular fins are 
assumed to have adiabatic tips. Results indicate the heat 
transfer rates for these two profiles vary by only 0.40 per 


cent for both tapered and cylindrical condenser. 
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ee Ne wm Ue LON 


pee THE ROTATING HEAT PIPE 

ite mOtacing, wiekless heat pipe 1S a closed container 
designed to transfer large amounts of heat from rotating 
Ma@eninery components. Essentially, it consists of three 
Main components: a cylindrical evaporator section, a con- 
Memser section which may be either tapered or cylindrical in 
Shape, and a fixed amount of working fluid. A typical tapered 
moeacing heat pipe 1S shown in Figure 1. 

When the heat pipe is rotated about its longitudinal axis 
mia Speed above a certain critical value, the working fluid 
forms an annulus in the evaporator section. Note in Figure l 
Meats the diameter of the evaporator 1s larger than the con- 
@emser. this larger diameter provides a greater liguid reser- 
tomer. AS heat 1s added to the evaporator, the fluid in the 
evaporator will vaporize. The vapor will flow axially towards 
me cOndenser as a result of a slight pressure difference, 
mealisporting the latent heat of vaporization with it. [In the 
condenser end, external cooling of the condenser causes the 
memonr tO condense. In the case of a tapered heat pipe, the 
centrifugal force due to the rotation of the pipe has a com- 
ponent acting along the condenser wall which accelerates the 


imequid condensate back to the evaporator to complete the cycle. 


Ie 





heat out 


Baer. lh. 





Schematic Drawing of a Tapered 
Condenser Rotating Heat Pipe 
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Percylindrical heat pipe, on the other hand relies on a hydro- 
mettic pressure gradient to drive the liquid condensate back 


to the evaporator. 


B. BACKGROUND 

ime stirst theoretical investigation into the performance 
of a tapered rotating heat pipe at the Naval Postgraduate 
School was accomplished by Ballback [Ref. 1] in 1969. He 
examined the limits in heat transfer controlled primarily by 
maid dynamic considerations. In particular, he considered 
the following four limits on heat pipe performance: a) the 
Bogling limit, b) the entrainment limit, c) the sonic limit 
aided) the condensing limit. Tantrakul [Ref. 2] calculated 
these limits for a specific heat pipe. He found the condens- 
mieerimit was the controlling limitation. In fact, the cal- 
culated heat transfer rate, based on the condensing limit was 
1/10th the heat transfer rate for the next lowest limit, the 
Seerainment limit. 

mieerdér tO Overcome thas condensing limitation and thus 
increase the heat transfer rate of the rotating heat pipe, 
mie concept of an internally finned tapered rotating heat 
Pipe was considered by Schafer [Ref. 3]. Schafer developed 
an analytical model for this tapered heat pipe with a trian- 
gular fin profile as shown in Figure 2. He assumed one dimen- 
eeonal heat conduction through the wall and fin. Corley [Ref. 


4] developed a two-dimensional heat conduction model using a 


ee 








Figure 2. Axially Finned Condenser Geometry Showing 
Fins, Troughs, and Lines of Symmetry 
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Finite Element Method formulation for the same geometry. To 
overcome the problem of few nodal points along the fin surface, 
Corley assumed a parabolic temperature profile along the fin 
emetace. JTantrakul [Ref. 2] modified Corley'’s computer code by 
increasing the number of finite elements from two to three in 
emder tO Minimize the error at the apex of the fin. Purnomo 
[Ref. 5] developed a two-dimensional Finite Element Method 
solution to the steady state heat conduction problem using a 
linear triangular finite element. Davis [Ref. 6] modified 
Purnomo's code to make it compatible with COPES/CONMIN [Ref. 8], 
an optimization program. Davis, in his modification, found a 
COding error in Purnomo's [Ref. 5] code, that once corrected, 
Bermitted Purnomo's corrected code to converge to Schafer's 
[Ref. 3] results. 

Purnomo's [Ref. 5] code is limited in that it is restricted 
imemone particular condenser geometric configuration, namely: 
aeeaxially finned tapered Condenser heat pipe with a triangular 
fin profile. In that tapered finned condenser are difficult 
to manufacture, it is doubtful that any widespread practical 
Mmeeeication of this geometric configuration will result. Cylin- 
drical condensers, on the other hand, can be manufactured with 
Meme less difficulty and might find practical application. 

This being the case, a more practical and beneficial code 
would be one that could analyze cylindrical condenser rotating 
heat pipes, both finned and smooth. In actuality, the most 


Deneficial code would be one that could analyze the following 
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mour geometric cCOnmfigurations: 1) tapered-internally finned, 
2) tapered-smooth, 3) cylindrical-internally finned, and 4) 

@ tindrical-smooth. The theoretical heat transfer performance 
Ot the four geometries could then be compared to determine the 
advantages and disadvantages of each design. An additional 
advantage would be gained if different fin profiles, i.e., 
merangular vice rectangular, could also be analyzed and 


eompared. 


meee labolS OBJECTIVES 
The objectives of this thesis are: 

ieee Develop analytical models for both cylindrical-smooth 
and cylindrical-axially finned condensers. 

2) Develop solution techniques to these analytical models 
that will account for temperature variations along the 
axial length of the condenser. 

3) Modify Purnomo's [Ref. 5] code to provide a solution to 
the two-dimensional steady state conduction heat transfer 
problem for the following four geometric configurations: 
a) tapered-smooth, b) tapered-finned, c) cylindrical- 
smooth, and d) cylindrical-finned. 

4) Modify Purnomo's [Ref. 5] code to provide the additional 
capability of analyzing a rectangular fin profile with 
an adiabatic tip. 

2) Obtain and compare results of the four geometric con- 


figurations given above for various operating conditions. 





ie TatGReELI GAL ANALYSIS FOR A CYLINDRICAL HEAT PIPE 


A. INTRODUCTION 

biea cyieacdricaieeondenser heat pipe, the radius of the 
condenser is constant along the axial length of the condenser. 
The flow of the condensate in the absence of vapor-liquid 
interfacial shear, is dependent upon the variation in hydro- 
static pressure with changes in film thickness along the sur- 
face of the heat pipe. Leppert and Nimmo [Refs. 8 and 9] 
mvestigated the phenomenon of film condensation on a flat 
horizontal plate. This situation is similar to film conden- 
meeron On the inside surface of a rotating cylindrical con- 
denser. In the case of a cylindrical condenser, the body 
merce, rather than being the force of gravity, is now the 
Mem@eritugal force caused by the rotation of the heat pipe. 
Weigenseil [Ref. 10] and Tantrakul [Ref. 2] compared experi- 
mental results for a cylindrical condenser rotating heat pipe 
With the theoretical results of Leppert and Nimmo [Refs. 8 
and 9] and found good agreement. The Leppert and Nimmo solu- 
tion was limited in that it was based on a constant surface 
mmperature along the length of the plate. A rotating heat 
pipe, in actuality, has a temperature variation along the 
axial length of the condenser which in some cases, may be 
Significant. This being the case, it was necessary to develop 


a mathematical model which would consider the axial temperature 





maeratron igetne Solution Of the heat transfer analysis. In 
the mathematical development that follows, a cylindrical 

smooth (unfinned) condenser will first be considered in that 
it is the simplest case. The model will then be extended to 


include a cylindrical axially finned condenser. 


B. THEORY FOR A CYLINDRICAL SMOOTH CONDENSER 
i Assumptions 
In developing the theoretical analysis, the following 
assumptions are made: 
a) Film condensation, not dropwise condensation occurs in 
the condenser. 
b) The condensate film undergoes laminar flow. 
c) Momentum changes through the condensate are small. 
Thus, there is essentially a static balance of forces. 
d) The vapor exerts no drag in the condensate; there is 
no interfacial shear. 
e) The temperature distribution within the film is linear. 
f) The vapor space is essentially at one pressure, oe 
g) The density of the fluid is much greater than the den- 
Sity of the vapor. Thus, the density of the vapor can 
be neglected. 
fee ihe centrifugal force is much greater than the force 
of gravity and, thus, gravity may be neglected. 
1) Velocity gradients in the circumferential direction 


relative to the pipe wall are negligible. 


Z2 





j) The condensate film thickness is much less than the 
EadiMiseOmeclUrVature Of the condenser wall. 

Ie Lhe rotating heat pipe is operating at steady state 
conditions. 


2. Condensate Momentum Equation (X-Direction) 


By applying the above assumptions and the coordinate 
eeeotem shown in Figure 3, an analysis similar to Nusselit's 
@riginal film condensation film theory may be used [Ref. 11]. 
Based on assumption c, a static force balance may be taken 


Sean infinitesimal fluid element in the x-direction as shown 


in Figure 3. This force balance results in the following 


equation: 
Die yO ee a6 (eqnez 2) 
X " Oy OX 
where +t = shear stress (l1bf/ft*) 
p = pressure (l1bf/ft* ) 
X = co-ordinate measuring distance along surface (ft). 
y = co-ordinate measuring distance normal to surface (ft). 


3. Condensate Momentum Equation (Y-Direction) 


Meomolitbaremaninen, Using Figure 5, a force balance 


Mmecne y-direction yields: 


LF =0: op wer = 0 (eqn 2.2) 
oy “i 


where P- = density of the fluid (lbm/ft?) 


W) angular velocity (rad/hr) 


i aeads (ft ) 





va 
oT 0 
Ty p+ 35 dy 





UTZ. Cr iT) 


Figure 3. Cross Section of Infinitesimal Fluid Element on 


Cylindrical Condenser Internal Surface 





eee acl WE te eaigny 


Integrating equation (2.2) between the limits y and 
é6* for y and corresponding limits of P and P| for pressure 


results in the following equation: 


Eee. OP pw°T (6*-y) (eqn 2.3) 
where ees Giesmmessumes ort the vapor (lbf/ft<), and 

hoe ee im tmackness (ft) 
Differentiating equation (2.3) with respect to x yields the 


following expression for dP/dx: 


dp _ 4p o-w*r dé* 
mee Ct aie (eqn 2.4) 
dx dx dx 


Applying assumption (f), Ue 1s constant, therefore, dP /dx=0) 


and substituting equation (2.4) into equation (2.1) yields: 


x 
ec p ,w°r = (eqn 2.5) 


V< 


Maeegrating equation (2.5) with the corresponding limits of 
integration y to 6* and t to 0 results in the following ex- 
pression for shear stress: 

2 ao: 


ob - §* R220 
fc dx Ly | (eq ) 


ZS 





t = —— (ean .1)) 


where vy = fluid dynamic viscosity (lbm/ft-hr) 


u condensate velocity (ft/hr) 
emibstituting equation (2.7) into equation (2.6) and integrating 
With the corresponding limits of integration 0 to y and 0 to u 


yields: 


[ i - yo*] (eqn 2.8) 


mile average velocity of the condensate may be found in the 


following manner: 














a Se 2 * 
l l P pw eGo 2 i 
i) = danke udy mae f i= y6 | ey (eqn 2.9) 
6 O 6 O u dx 
or 
e. Pew°r dé* 5. 
0 AE [+] (eqn 2.10) 


Eee CONtIinuity Equation 


The continuity equation for mass flow requires that: 


ieee uA bean 2.11) 


where m condensate mass flow rate (lbm/hr) 


A = cross sectional area of the fluid (ft?) 
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fii1s can also be written as 


* 


6 - 
m -{ Pu 2urdy (eqn 2.12) 
O 


Substituting equation (2.10) into equation (2.12) and inte- 
grating yields: 


Zi 


Zee * x3 
ener Olle Ore a = (eqn 2.13) 
u 


Differentiating this equation with respect to x yields: 


: a ,d6* 6%, 
dx u dx dx 5 





(eqn 2.14) 


O ISIMe ms oy: Equation 


Having applied assumption (e), if the film surface 


temperature is at the saturation temperature (T ) of the 


sat 
vapor and if the wall of the axial increment 1s at a given 
constant temperature Glee) then the heat transfer by conduc- 


mom Of a fluid element of surface area dA is: 


ke (Tsat - Tw) dA 
eae 


aah. = F (fecie 25. leap) 
Where dq = differential heat transfer rate (Btu/hr) 

lame ircdoctit- ) 

ke = thermal conductivity of the condensate film 


Cie Neat te Pp ) 


a 





ie 


ave saturation temperature (degrees F) 


T,, = inside condenser wall temperature (degrees F) 
Considering the change of phase and defining he, as 

© 
the average enthalpy change of the vapor in condensing to a 


liquid and subcooling to the average liquid temperature of 


the film, then dq is also defined by: 


dq = eg din (eqn 2.16) 
where Neo = laten heat of vaporization (Btu/1bm) 
Cy = specific heat (Btu/lbm R) 
AT = Cla - ee) 
Ore = Nee tamOr, 55. 7° Cy sae 


Rearranging equation (2.16) and substituting this equation into 


equation (2.15) yields: 


dm . Ke(Toat - Ty) 27r 
ee 


ae 5 feqn™2e7,) 


Finally coupling the energy and continuity equations result in 


the following differential equation: 


x dé* *& 3 < = ’ 
6 — - 6 j]- —— Ty) - (eqn 2.18) 
eG 


Equation (2.18) can be solved using the Finite Element 


Method to provide the film thickness profile along the axial 


28 





Memeth of a cylindrical condenser. Appendix A provides a 
ieeriied description Of this solution. Once the film pro- 
file is known, a steady state two-dimensional heat conduc- 
tion analysis can be performed. 
7. Determination of Heat Transfer Rate 

Assume that the cylindrical condenser section of the 
rotating heat pipe is divided axially into a number of incre- 
ments. Then for any axial increment of a cylindrical conden- 
ser, the differential heat flux can be determined by the 


molilowing expression. 


dar = eS (eqn Dos) 
6 4 JUL ae 
Be : 
W ext 
where To = ambient temperature (degrees F) 

thick = thickness of the condenser wall (ft) 

Se = thermal conductivity of the wall material 
Chey Mae t= ) 

ieee = external heat transfer coefficient 


(Btu/hr-ft?-F) 
Note the three terms in the denominator are the thermal re- 
Sistances of the film, wall and external convection respectively. 
The differential heat transfer rate for any increment 


can be found by the following relationships: 


dqe= dq)’: 2nrdx (eqn 2.20) 






Or 

Pees iat ae we x 

a Sa ea 
me Geka 1 


fF W ext 


(eqn 2.21) 


Equation (2.21) represents the total heat transfer rate for 
Bieincremental section of width dx. To find the total heat 
transfer rate for the entire cylindrical condenser, the in- 
cremental heat rates must be summed over the entire length 
of the condenser. Therefore: 

NDIV 


= ¢£ dq Ceana2. 22) 


Q 
Oecd 7 


where NDIV = total number of axial increments. 


G. THEORY FOR A CYLINDRICAL AXIALLY FINNED CONDENSER 
mo ASsumptilons 
Referring to Figure 4, it is obvious that the analysis 

eeeaecylindrical internally finned condenser is more compli- 
cated due to the mass flow from the fins into the trough 
region between the fins. For this reason, in addition to 
the simplifying assumptions made for the smooth condenser 
which are listed in the previous section, the following 
assumptions must also be made: 

a) Referring to Figure 5, the mass flow along the fin sur- 


Picee@ocs MOL tlOw axidadkty in the x-direction, but only 





D) 


€) 


Pom EE Gemourrace. Or che fin in the Z-direction into the 
POM inNUSeiidsserlow In the axial direction is only 
permitted in the trough region between the fins. This 
1s a reasonable assumption in that the film thickness 
along the fin surface is very small in relation to the 
film thickness in the trough. This being the case, the 
MdrOStatiG temec in the x-direction on a fin fluid ele- 
ment will be much less than the centrifugal force com- 
pementeim the Z2=-direction om that same fluid element 
Pomeciigmtiat t£lutd element into the trough. 

Suctmacwim the axial direction, there 1S no pressure 
eEMancemalona the Surface Of the fin in the z-direction. 
E_ieeepeeassumed that the temperature along the’ con- 
vective surface of the fin is at a constant value ave 
This average fin surface temperature is the arithmetic 
ayenragesOf the fim tip temperature and the fin base sur- 
face temperature where the fin intersects with the wall 
of the condenser. This is a valid assumption if the fin 
section is divided into a sufficient number of finite 
elements. Purnomo's [Ref. 5] results indicate a less 
than one degree variation, even for very large fin half 
angles. This variation in temperature will have an in- 
Significant effect on film thickness along the surface 
of the fin and can be neglected by using an average 


value. 


on 





Figure 4. 





Trough Section of Axially Finned Condenser 
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Figure 5. Axially Finned Condenser Symmetric Fin Section 
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a, UMass Plow in the X-Direction 
As a result of assumption (a), the resulting momentum 
equation in both the x and y directions as well as the equa- 
tions for velocity and mean velocity are identical to those 
developed in section B for the smooth condenser and will not 
be redeveloped here. Looking now at mass flow in the x-direc- 
tion which is limited to flow only in the trough, the mass 


tlow rate 1S given by the following expression: 


2,.2 - 
eer = — g*°(c6* + 5** tana) (eqn 2.23) 
u x 


i) 


where a fin half angle (radians) 


I 


E Wid cime te tne urougn (tt) 

Meee that the quantity in parentheses is the cross sectional 
area of the film condensate in the trough (See Figure 4). 
meerme the derivative of equation (2.23) with respect to x 


fMeetas the rate of change of mass flow in the trough for a 


given axial increment. 


° 2 
dm. otal a a d ‘e 
7 lax Boot 
Ging Su dx . 


+ §**tana) ] (eqn 2.24} 


Equation (2.24) represents the rate of change of the 
total mass flow rate with respect to x in the x-direction. 
This equation must be coupled with the energy equations for 
the fin and trough to develop a representation of the film 


Peorile in the trough. 





S.. 


Massetlow im the Z2Z-Direction 


Examining an infinitesimal fluid element on the sur- 


face of a fin for any axial increment of width Ax, as shown 


in Figure 6, the momentum equation in the z-direction becomes: 


where 





OT 

al _3P 2 
sy : P pw* TCOSa (eqn 2.25) 
Peeemotcarestuess INtne 2-direction ios a iecea) 

z = co-ordinate measuring distance along the surface 


of the fin (ft) 


Neglecting dP/dz based on assumption (b), and integrating 


equation (2.25) from T 


where w 
6 
Note, 6 


me. the 


, to 0 and y to 6 yields: 


T = ie = p ew rcosa(s-y) (eane 2.20) 





willbe velocity Im the z-direction (ft/hr) 


fin film thickness along the surface of the fin (ft) 
, the fin film thickness should not be confused with 


film thickness in the trough. Integrating equation 


(2.26) from 0 to w and 0 to y provides the following expression 


meet iuid velocity: 


9-6" rcosa 


y? 
Ww = =e (d(z)y - oe) (eqn 2.27) 


Se) 








dy 
on 
Z U/ at sy dy 
OL 
dz 
Figure 6. Cross Section of Infinitesimal Fluid Element on 


Fin Surface of Axially Finned Condenser 


He 





This relationship may be used to find the average fin fluid 


velocity w: 


5 0 Ww rcosa 6 (2 )i- 
i = J wdy = ————————_ (eqn 226) 
OQ 


SU 


Or fe 


Mie mass flow rate in the z-direction along the surface of the 


fin for a given axial increment is given by: 


ee P -wdA (eqn, (2.29) 


where dA = the cross sectional area of the fluid flowing along 
the fin surface (ft*) 
eeestituting equation (2.28) into equation (2.29) yields: 
0 - w° Tcosa §(z) dx 


= a eee Geant2 450) 


m.-- 
forn a 


This equation is identical to equation (15) of Schafer's [Ref. 
3] analysis if the condenser cone half angle (ff) is set equal 
to 0 which is the case for a cylindrical condenser. 
ee Energy Equation for the Trough Condensate 
An energy balance on an infinitesimal fluid element in 
the trough of an axial increment of width dx with surface area 
€ - dx yields the following expression for heat transfer by 


Sonauction: 


eee T - T )edx 
44+ rough fe osat W 
6 


(eqns 2 = 51) 
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Mere also, that the trough heat transfer rate is given by: 


dd + rough 7 ae u 


M+ rough (ean 2452) 


Sembining equations (2.31) and (2.32) and dividing by dx re- 


sults in the following: 


dm rte cea) 
EEoueh | & a W amen) 
dx Reo 6 


Bauation (2.33) is an expression for incremental change in mass 
flow rate with respect to x due to condensation in the trough 
region. 
a. Energy Equation for the Fin Condensate 
An energy balance on a differential element of surface 
area dx-dz yields the following relationship for differential 


heat into the fin: 


Nee i eee 2) ax 2 
fees a te i eT 
fig. = h di-- = ———__ (eqn 2.34) 
fin fg fin 5(z) 
where ean (z) = fin surface temperature at some position z 


alemomete Suriace Of the fin (degrees F) 
Since the fin condensate mass is assumed to flow only in the 
@eetrection, equation (2.30) is differentiated with respect 
Bemzeand Substituted into equation (2.34). After substitution 


and rearrangement, the following equation results: 
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k-(T... - Te:.(z))dz 
$(z)3dé(z) = = _s8t tin (eqn 2.35) 
oes Neg cosa 


AMivinicmasisumptilon (Cc) 1.e., T e367) equals lave one 


miele atid integrating equation (2.35) from 0 to & and 0 to z 


frelas the following relationship for fin film thickness 4 (z): 


1/4 
4 k.(T - T aul 
§(z) = oe (eqn 2.36) 
OP, wT Neg cosa 
where ee = average fin surface temperature (degrees F). 


Suestituting equation (2.36) into equation (2.30) and solving 
for rate of change of mass flow rate of the fin with respect 


woex £LOr an increment of width dx yields: 


1/4 
2 Do 4 a % 
dm; _ 20 ¢ w°rcosa - aerate nasi Le amen 
eee 
ax oul OP, wT Nes cosa 
Mere z* = z-d6*/cos(a). Note, z* is the distance along the 


Mumeace Of the fin from the fin tip to the trough film thick- 
ness (5 *). Note also, that the right hand side of equation 
mo’) 1s multiplied by two; this accounts for mass flow from 
Mieertins on both sides of the trough. 
eee Continuity Equation 
For any axial increment of length dx, continuity dic- 


mates that: 


oS. 





dm Giiileae dm 
Ota! oan trough - 
cao oumace dx ceo) 


m@bstituting equations (2.24), (2.335) and (2.3/7) into equation 


(2.38) and rearranging yields: 


Scar Lie 
meee 2S” Ce gkts §**tand)] = - —t Sat Ww 
Boo dx o-2w?r A 
if fg 
; . 3/y 
keh - T NZ 
- 25*cosa | i s8t sve (eqn 2.39) 


oe 
oF Sc Neo cosa 


memaeton (2.39) can be solved using the Finite Element Method 
Pormulation provided in Appendix A. The solution of this equa- 
tion provides the film thickness profile along the axial length 
of a cylindrical finned condenser. 
7. Determination of the Heat Transfer Rate 
Once the film profile has been determined within the 
trough, the local convective heat transfer coefficient can be 


found for the trough using the following relationship: 


“s 


h(x) a 
trough ale 


(eqn 2.40) 


In a similar manner, the local heat transfer coefficient along 


the surface of the fin can be found by: 


k 
1 
A(z) 3, - Sz) (eqn 2.41) 
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fieeditterential heat transfer rate for any fin section, as 
shown in Figure 4, of axial incremental length dx is: 
metor the trough: 


; Chae - TJ €dx 


cele ough: ~ ao... (eqn 2.42) 


rough 


where €°dx is the surface area of the trough, and 


be tor the fin surface: 


Zo. «((T =k )dxdz 

ee (eqn 2.43) 

aT 
z NZ) ean 

where zis the surface length of the fin. 

Mies total differential heat transfer rate per axial increment 


meecound by summing equation (2.42) and (2.43) for the total 


mumper ©: fins. That is: 
NFIN 
d4-otal 7 : (ddein : 14+ rough? (eqn 2.44) 


mavere NFIN is the total number of axial fins. 
faea Similar manner, the total heat transfer rate for the entire 


finned condenser can be found by the following relationship: 


NDIV 
Ore 7 , ae oneal (eqn 2.45) 


where NDIV is the total number of axial increments. 
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A. 


Ci COMPU NeCODE me DEoGR LPT FON 


GENERAL DESCRIPTION OF CODE 


The computer code consists of a main body and eight sub- 


routines. Basically, the code which is provided in Appendix 


Peis @ modification of Purnomo's [Ref. 5] code. The function 


of each subroutine used in the code is as follows: 


a) 


Db) 


@) 


d) 


"CORRES" established the correspondence between the 
local and global nodal points used in the Finite Element 
Method solution for the two-dimensional steady state 
Meiumconanculoneproentem. In so doing, "CORRES" also 
numbers all elements and nodal points inthe finite ele- 
ment model and assigns local nodal points to each of the 
Pecicnec wee nedddition “CORRES” also defines major ele- 
ment numbers used in other subroutines as control 
parameters. 

"COORD" defines the x and y coordinates for all nodal 
points in the finite element heat conduction problem 
model. 

“"DLSTAR" determines the film thickness (6*) on the sur- 
face of a smooth condenser or in the trough in a finned 
condenser. 

PinGokn “determines the heat transfer coefficient for 


all convective surface elements. 
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e) ''FORMAF'' formulates the Finite Element Method equations 


for the two-dimensional steady state heat conduction 


ao bake. 
im BANDEC' is an equation solver for a symmetric matrix 
which has been transformed into banded form. ‘'BANDEC" 


will return the solution to the two-dimensional heat 
conduction problem. 

g)- “HTCALC" determines the elemental, incremental and total 
heat transfer rates. 

h) '"DELCRV'" determines the condensate film profile in a cyl- 
indrical condenser. 

Two additional Naval Postgraduate School computer library 
routines are also used in the code: 

a) "DPOLRT'" is a nonIMSL double precision library routine 
that determines the roots of a real polynomial. This 
routine is called by "DLSTAR" to determine the film 
thickness for the succeeding increment in the analysis 
of a tapered condenser. 

oe LEQTZE" is an IMSL double precision libary routine that 
sOlves a set of simultaneous linear equations. This 
routine is called by “DELCRV" to solve the Finite Element 
Method equations for the cylindrical condenser film pro- 
file problem. The resulting film profile is then used 
in the heat conduction analysis. 

In order to use the computer code to analyze heat transfer 


in a rotating heat pipe, nine data cards are required. A user's 





guide describing these data cards and required input is pro- 
vided in Appendix B. The input data, describing the geometric 
Seamriguration of the rotating heat pipe as well as the operating 
parameters determines which solution technique is utilized in the 
analysis. The solution technique for each of the four con- 
denser geometries, l1.e., tapered-smooth, tapered-axially finned, 
mamindrical-smooth and cylindrical-axially finned is different. 
In all cases however, the Finite Element Method is used to 

solve the two-dimensional steady state heat conduction problem. 
This solution is the one developed by Purnomo [Ref. 5] and 

has not been modified. Details of the development of this 
solution are described in detail in Purnomo's thesis [Ref. 5] 
and will not be repeated here. This being the case, each of 


the four solution techniques will now be discussed in detail. 


ieee NTERNALLY FINNED TAPERED CONDENSER SOLUTION 

The complete development of this solution technique is 
provided in Purnomo's [Ref. 5] thesis and will not be rede- 
veloped. When an equation is required for clarity, the equa- 
tion in final form will be provided. Where there is a 
ieart ication to Purnomo's [Ref. 5] code, this modification 
Will be noted. 

The condenser of the internally finned-tapered condenser 
1s divided into NDIV axial increments. These axial incre- 
Memes are then suddivided circumferentially into ZFIN number 


of subincrements where ZFIN is the total number of fins. 
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These subincrements are then divided in half to form the 
mrsic Symmetric unit, one-half of a fin-trough section as 
shown in Figure 1. This unit section is then divided into 
a number of linear triangular elements. The number of ele- 
ments depend on the input parameters. The only limitation 
is that the number of system nodal points must not exceed 
109; otherwise, certain variables, e.g., x and y, would 
exceed allotted storage values. Figure 7 shows a unit 
section subdivided into 25 elements. After this unit is 
subdivided, each nodal point 1s assigned an x and y coordi- 
nated based on the geometric input parameters. 

To start the iteration, two initial values are required: 
1) an initial temperature for the nodal points along the 
internal convective boundary, and 2) an initial trough film 
thickness (6*). The initial temperature is provided as an 
input parameter and the initial trough film thickness at the 
first increment is provided by a relationship taken from an 
analysis by Sparrow and Gregg for condensation on a rotating 
@ick [Ref. 12]. 

Once these values are known, the heat transfer coeffi- 
Clent for the internal convective surface elements are found 


using the following relationships: 


1/4 


). = | ==—_,—___58 (eqn 3.1) 
in ea) 4u(-AA+z /3-BB+2°/2+(T.,,-T,)2 


Kk. Keo -°w* +xsinf he cosfcosa 
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Fisure 7. Axially Finned Condenser Symmetric Section Sub- 
teeidedmiimtomsoeLinear Itiangular Finite Elements 





where x = distance from condenser end (x=0.0) to midpoint 
of increment (ft) 
i= fin apex temperature (degrees F) 
AA, BB = constants in the parabolic temperature determined 


by a Langrangian fit. 
memeche trough surface elements: 
ae 


hOD trough = S¥(x) (eqn 3.2) 


These heat transfer coefficients, along with the thermal 
conductivity and x and y coordinates of the nodal points are 
used to form the Finite Element Method equations for the two- 
dimensional heat conduction problem. The equations are then 
sOlved to yield a temperature distribution in the symmetric 
Secrtion. 

The above iteration is repeated, where, now the solution 
memperature distribution from the previous iteration is used 
Memcalculate the heat transfer coefficients along the con- 
vective surface fin elements as well as a new 6* using Sparrow's 
mmaearege’sS relationship [Ref. 12]. This new 6*, in turn, is 
Meeq tO determine the heat transfer coefficients of the trough 
Eements. 

Again, the Finite Element Method subroutines will yield 
meccmperature distribution for the symmetric section. At 
this point, the nodal point temperatures are checked for con- 


vergence using the following relationship: 
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Meee |) 2 Loj-1 | < crRIT i= 1,2,...NSNP (eqn 3.3) 


where NSNP = number of system nodal points 


j = present iteration 
j-1 = previous iteration 
CRIT = convergence criterion 


If this convergence test is successful for all nodal points, 
the increment is considered solved. If any one nodal point 
fails, the iterative process is repeated until convergence is 
met. The convergence test in equation (3.3) is different than 
the one used by Purnomo [Ref. 5] in his thesis. Purnomo com- 
feeeed incremental heat transfer rates per unit of condenser 
length, Q;> rather than temperature as is done in the modified 
code. 

If convergence is met, the heat transfer rate is determined. 
From this heat transfer rate, the incremental mass flow rate is 
determined by the following equation: 


2Q5Ax 


Me otal ~ he (eqn 3.4) 
Oo 


where Q, = heat transfer rate per unit length (Btu/hr-ft) 


Ax incremental width (ft) 
Using this value of incremental mass flow rate determined 
by equation (3.4), the following equation is used to calculate 


the subsequent interval's trough condensate film thickness (§*) 


With a polynomial rootfinder subroutine: 
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. Pw” (r+xsinf) 6*? (x)sinp(6* (x) e+6** (x) tana) 
meotal « 3y 





(eqn 3.5) 


where € = trough width (ft) 

This resulting value of 6*(x) is then defined as the trough 
film thickness for the next increment. In addition, the solu- 
mmom temperature distribution from the previous iteration is 
used as the starting temperature distribution for the next 
increment. 

its iterative process at each increment is repeated until 
convergence is met, and, is continued at each increment until 
j@@ewentire length of the condenser has been transversed. In- 
cremental heat rates are then summed to yield the total heat 


Meanster rate. That is: 


= JR? * ° 3 
otal 2°ZFIN* 2 Q, Ax (eqn 3.6) 
where ZFIN is the total number of axial fins. 
Once the total heat transfer rate has been determined, the 


meoolem 1s solved and pertinent data is provided as output. 


C. SMOOTH TAPERED CONDENSER SOLUTION 

The heat pipe condenser is divided into NDIV number of 
axial increments as in the finned-tapered condenser solution. 
These axial increments are then subdivided into 360 segments 


of equal length; these segments are the basic symmetric unit. 
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mrss unit section. 15 divided into a number of linear tri- 
angular elements with the same limitation as before; the 
number of system nodal points must not exceed 100. The 
system nodal points are then assigned x and y coordinates 
based on the input geometric parameters. 

Wiomstart the iterative process, as in the finned-tapered 
case, the initial value of temperature which 1s an input 
parameter is used to solve for the initial value of fin 
thickness (6*) based on the Sparrow and Gregg analysis [Ref. 
|. 

Omee fhis imitial value of 6* is known, the heat transfer 
coefficients for the internal convective elements can be 
determined using equation (3.2). These heat transfer coef- 
[Meerents are used in the Finite Element Method equations. 

Mie E€quations are solved yielding a temperature distribution. 
Piesiteration is repeated until convergence is met, just as 
ietene finned-tapered case. 

When convergence is met, that is equation (3.3) has been 
Satisfied, a new film thickness §*(x) can be found by one of 


the following equations for §*(x): 


ie eee a. ue r z 
ete | > | = Dr, | cum cz - 1 = 0 Ceqmro. 7) 
or 
sta 
Sh Te a= 0 (eqn 3.8) 
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‘em eee 


Bketclaee- i) “/3 
Xx) ee ee [1 - {+} ] (eqn 3.9) 


20 -wrsin Phe, (r+xsinf) 


il 


O° (w*r-g)sinf hex” 


where Sh 
sat 


O-Tt he x* cos 
I ig ac 


Uke (Toot ; Ty? 
P eV cos@ 
Rev = 
u 
g —eeceteration due to gravity (ft/hr-) 
V = vapor velocity (ft/hr) 
ee = local wall temperature (degrees FO 
= = shear stress vapor-liquid interface (l1bf/ft?)} 


meumation (3.7) defines the film thickness distribution for 
Meemooth tapered rotating heat pipe derived by Daniels and 
feeeJumaily [Ref. 13]. This equation takes into account the 
eae effects of counter-flowing vapor. Equation (3.8) is a 
modification of Equation (3.7) neglecting the drag losses. 
Equation (3.9) was developed by Ballback [Ref. 1]. This equa- 
meen also neglects drag. 

Depending on a particular control parameter which is part 
Of the input data (See Appendix B) one of these equations is 


Heed to solve for the film thickness (6*(x)) for the next 
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Gncrement. In addition, the solution temperature distribution 
from the previous iteration is used as the starting temperature 
Seer i bution for the next increment. 

This iterative process at each increment is repeated until 
convergence is met, and is continued at each increment until 
the entire length of the condenser has been transversed, just 
mein the finned-tapered case. Total heat transfer rates are 
then determined by summing the incremental heat transfer rates 
for the entire length of the condenser by the following 
relationship: 

NDIV 


= 360* f Q; Ax (eames 6) 


Sora eal 


feo MOOTH CYLINDRICAL CONDENSER SOLUTION 

As in the smooth-tapered case, the condenser is first 
divided axially, then it is divided circumferentially into 
560 segments of equal length. These segments are the basic 
Symmetric unit section to be considered. Again the symmetric 
unit is subdivided into linear triangular elements and x and 
y coordinates are assigned to the system nodal points. 

To begin the iteration, an initial temperature estimate 


Which is an input parameter is assigned to the convective sur- 


face nodal points. Using this initial temperature estimate, 
the maximum film thickness, ae which is located at x = 0, 


j determined. This maximum film thickness value is then 


5/4 





used as one of the boundary conditions in the solution of 
@emetion (2.18) using the Finite Element Method. Equation 


mes) iS repeated here for reference. 


* d * *% 3 


5 Z fe 8) - _ ke (Tsat - Twit (eqn 3.11) 


Pp wr Aggy 


[mentinite element solution of Equation (3.11) provides the 
film thickness profile (6*(x)) at the midpoint of each incre- 
meme along the length of the condenser, but only the first 
mierement value is applied at this point of the analysis. 

Once 6*(x) is known at the first increment the heat transfer 
@eerticients for the internal convective surface elements can 
De determined using equation (3.2). The steady state heat con- 
duction problem is then solved and a temperature distribution 
Mpecne Unit section results. This process is now repeated. 

A new maximum film thickness, film profile and thus §*(x) at 
the first increment is found based on the solution temperature 
G@@stribution from the first iteration. With this new value for 
6é*(x), the heat conduction problem is again solved for a new 
temperature distribution. This iterative process, involving 
Seemciton of the film profile with each iteration is continued 
until temperature convergence is met at the first increment. 
When convergence is met, the film profile determined on the 
iteration in which convergence was met is used to provide the 
Maes Of 6*(x) for the remaining increments along the length 


Of the condenser. 
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Using this predetermined value of 6*(x), the iterative pro- 
cess 15 continued at each increment until temperature conver- 
gence is reached. As convergence is reached at each increment, 
the internal wall temperature, which is the same for all inter- 
mal wall nodal points of the section, is stored for future 
moplication. 

When convergence is reached at the final increment, equation 
(3.11) is once again solved for the film profile. It should be 
noted, however, that the right hand side of equation (3.11) is 
temperature dependent. It should also be noted that the temper- 
ature varies axially along the length of the condenser. In 
meager tO account for this temperature dependence and temperature 
variation, the right side of equation (3.11) is now determined 
for each increment using the wall temperatures that were stored 
Mmmmeacn increment. The finite element solution of the film pro- 
Tile will now account for the temperature variation along the 
bength of the condenser. 

This final film profile provides the value of §*(x) for each 
increment. The iterative process of solving for heat transfer 
coefficient, temperature distribution, and temperature conver- 
gence is continued for each increment until the final increment 
is reached. 

Bauation (2.13) provides a relationship for total mass flow 
rate as a function of position. The total mass flow rate at 


the overfall into the evaporator is given by: 
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2 2 * 3 


* 
nm = - anee ees == : St) cSatimo. 12) 


where §*(L) and dé*(L)/dx are the values for the film thickness 
meeerate of change of film thickness with respect to x at the 
evaporator end of the condenser. Another relationship for mass 
flow rate based on the steady state heat conduction solution 


meeegiven by: 


NDIV 
m= 360% 5 Sa0X (eqn 3.13) 
je) Neg 


im, in fact, the solution of equations (3.12) and (3.13) are 
equal, then the mass flow rate of the condensate returning to 
the evaporator is equal to the mass flow rate of the vapor 
being condensed on the surface of the condenser. Or, to put 
lt another way, continuity is satisfied. 

It should be noted that the film profile maintains the same 
Dasic shape, that is, 6¥*(x) at x=0 is always equal to Oe and 
decreases to a specific minimum value at x=L. This being the 
case, if the maximum film thickness is varied, the film thick- 
nexx profile will vary in the same manner. For example, if 
the maximum film thickness is increased, the entire profile 
Will also increase. This will result in an increased internal 


thermal resistance and thus lower heat transfer rate. As a 


result of the lower heat transfer rate, the mass flow rate as 
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determined by equation (3.13) will be less. At the same time, 
Mewever, the greater film profile will result in a greater 
femme Of S6“(x) at the overfall, x = L. Yet, since the profile 
maintained the same basic shape, the derivative at the overfall 
remains relatively constant. Thus the mass flow rate as deter- 
fmemed by equation (3.12) will increase. A decrease in the 
feemmum film thickness will result in an Opposite effect to 
the mass flow rates determined by equations (3.12) and (3.13). 
This being the case, the mass flow rates, as determined 
by equations (3.12) and (3.13) are now compared to determine 
feeene film profile is in fact the solution profile to the 
problem. If continuity is not satisfied, the maximum film thick- 
ness is varied and the entire iterative process is restarted. 
This process is continued until the film profile mass flow rate, 
equation (5.12) converges towards the heat transfer mass flow 
rate, equation (3.13). When the absolute difference between 
these mass flow rates is less than a specific value, the re- 
Sulting heat transfer rate is considered the solution to the 


problem. 


ier i NNED CYLINDRICAL CONDENSED SOLUTION 

As in the finned-tapered case, the condenser is divided 
axially and circumferentially into the basic symmetric section 
as shown in Figure 1. This unit is then subdivided into linear 
triangular elements and x and y coordinates are assigned to 


Seen nodal point. 


56 





homie ctmerne ILeraulye mnocess, as before, an initial temp- 
erature estimate is assigned to each nodal point along the in- 
ermal convective surface. An initial trough film profile is 
determined, using equation (3.11) and this initial temperature 
estimate. In this case, however, the maximum film thickness 


es* 


a) is not calculated but is an input parameter. 


Once §6*(x) is known at the first increment, the internal 
heat transfer coefficients are determined, using equations (5.1) 
and (3.2). These values are then used in the Finite Element 
Method solution of the steady state heat conduction problem. 
Becemperature distribution is determined and the iteration is 
repeated until temperature convergence is met. Note that a 
Mew film profile is determined for each iteration. 

As in the cylindrical-smooth condenser case, once conver- 
gence 1S met at the first increment, the film profile that was 
determined for the iteration prior to convergence at the first 
mierement is then used to provide the film thickness §*(x) for 
the remaining increments. 

merthe final increment, a new film profile for a finned 
cylindrical condenser is then determined by solving the follow- 


ing equation developed in Chapter II: 





* 
- ~ ss (g* e+ §**tano)] = - Sk ¢(T t_- Tw) ue 
Oecu-r i 
4k .(T - T Nae: 3/4 
= 26*cosa a once) 


PO, wr he, cosa 
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A description of the solution of this differential equation 
Using the Finite Element Method is provided in Appenxix A. 
The iterative process of finding the solution temperature 
@stribution for all increments is then repeated until the 
length of the condenser has been transversed. 
[ethics point, the total mass flow rate of the condenser 


returning to the evaporator is determined by the following 


relationship: 
Op wer dsg™ 3 4 
m air ae Cag 45 tank ) ZFIN (eqn 3.15) 


where ZFIN is the number of axial fins. 
ios mass flow rate 1s compared to the mass flow rate given 
by the following equation: 

NDIV 


wh = ZFIN*2* & Roe (eqn 3.16) 


i=l ee 





HNst as in the smooth-cylindrical condenser case, if the 
absolute difference between the two mass flow rates is less 
than a mass flow convergence criterion, the problem is con- 
Sidered solved. If not, Sa 1S varied and the entire itera- 
tive process is started again. As nthe smooth-cylindrical 
Semaenser, varying ors Will have the same effect on the film 


profile and the heat transfer rate. Since the temperature dis- 


tribution along the condenser has been solved once and closely 
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approximates the final solution to the problem, on successive 
iterations, equation (3.14) may be used to solve for the film 
mgerile rather than equation (3.11). Equation (3.11) was used 
on the first iteration because the finite element solution con- 
verges more quickly than equation (3.14) when an estimated 
temperature is used. 

A word of caution 1s required. The solution of equation 
(3.14) is highly sensitive to the value of a ht ethe™ initial 
value of Sate is inconsistent with the actual solution, e.g. 
too small, the Finite Element Method solution of the film pro- 
Mmee will not converge. If this 1s in fact the case, the prob- 
lem will be automatically terminated. A new value of OF as should 
then be chosen and the problem restarted. | 

One additional topic which should be addressed is the rec- 
tangular fin solution process. The rectangular fin profile is 
a Slight modification of the finned-tapered or finned-cylindrical 
Condenser solution. The only variation is that the top elements 
meeetie rectangular fin are assigned heat transfer coefficients 
Sees Oo. thus, the tip of the rectangular fin is considered 
adiabatic and no heat is transferred through this top face. 
Other than this modification, the solution techniques are the 


Same as tor the finned cases addressed above. 
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iv> RESULTS AND DISCUSSION 


Prior to the heat transfer analysis of the various con- 
wemser cOnfigurations, it was mecessary to verify the finite 
element solution of the film profile. The development of this 
Somution is discussed in Appendix A. 

Leppert and Nimmo [Refs. 8 and 9] had developed an analy- 
fieal sOlution for film condensation on a horizontal plate at 
a constant surface temperature. Their analysis and resulting 
differential equation is identical to the development in Chapter 
{Il for a smooth cylindrical condenser if the acceleration due 
tO gravity is replaced by a radial acceleration term. This 
being the case, this modification was made and the analytical 
solution was used as a reference for comparison. 

For the test runs of the finite element solution, a constant 
Surface temperature was assumed. In addition, a value for the 
maximum film thickness Cone and minimum film thickness (875 ) 


at the overfall were required. The value for ae was determined 


re 


based on a relationship developed by Leppert and Nimmo [Ref. 8]. 


The value for $F was arbitrarily chosen. 


Nh 


For identical geometry, surface temperature and maximum and 


minimum film thicknesses, the results of both analyses were 


identical. In order to develop confidence in the finite element 
2 x : * * 
Sevution, ein was Mamted trom 0,10 - ee Orta s a> See The 


resulting profiles agreed at all locations along the length of 
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the condenser. However, when the temperature was permitted to 
meay along the length of the condenser in the finite element 
solution and a corresponding average temperature was used in the 
maalytical solution of Leppert and Nimmo [Refs. 8 and 9], the 
profiles were no longer in agreement. This was to be expected, 
Meeericularily in the case where there was a sixteen degree 
Fahrenheit variation along the length of the condenser. This 
substantial temperature variation resulted in a significant 
fmeration of fluid properties which would account for a differ- 
ence in film profile. 

Due to the agreement between the finite element solution 
and the analytical solution of Leppert and Nimmo [Refs. 8 and 
9] for a constant surface temperature, it was decided that the 
mate element solution does provide a satisfactory representa- 
memo r the film profile. This being the case, the finite 
element solution was then incorporated into the code to provide 
the film profile for the cylindrical condenser. 

See the finite element solution of the film profile was 
verified, the heat transfer analysis could be accomplished. 

The analysis considered both copper and stainless steel con- 
densers with the following four geometries: a) tapered-smooth, 
b) tapered-axially finned, c) cylindrical-smooth and d) cyl- 
indrical-axially finned. Table I lists the geometric para- 
meters held constant for all analyses. In all cases, the 


working fluid was water. 
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TABLE I 


Condenser Geometric Parameters 
Held Constant During All Analyses 


condenser length 8.500 ineltes 


minimum radius io S75 inches 


Welechilekness On O Sins Pnehies 

inpeaddition, the following geometric parameters were also 
utilized when required. This requirement was based on the 
condenser geometry being considered, l.e., tapered-axially 


finned. 


TABLE II 


Condenser Geometric Parameters 
Applied as Required 


height of fin Or Grs20'o inches 


fin half angle = ZOOS degrees 
condenser cone 
half angle = 1.00 degrees 


In the analysis, the heat transfer rate was determined for 
the four different geometries listed above for both copper and 
Stainless steel. The ambient temperature was set at UOhe: 
and the heat transfer rate was determined for each possible 


combination of the operating parameters given in Table III. 
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TABLE ITI 


Operating Parameter Matrix 


Rotational Heat Transfer Dee ugelieskoyal 
speed Coecruverent Temperature 
(RPM) (pu tir-£t--F) (degree F) 
70 Ure VOOe 90.0 

1400.0 5000 LZ ORO 
2800.0 FO CO 150.0 
SO. 0 


Thus, for each condenser geometry, there was a total of 72 
mmaiyses, 36 for the case of the condenser with a copper wall 
and 36 for the case of the condenser with a stainless steel 
wall. 

The first condenser geometry considered was a smooth con- 
memoer. Figures 8-13 compare the heat transfer rates of smooth 
merindrical condensers with those of smooth tapered condenser. 
meparticular, Figures 8,9, and 10 indicate the results of the 
mmelyses Of smooth copper condensers at rotational speeds of 
700, 1400 and 2800 revolutions per minute(RPM) respectively. 
Meeures 11, 12, and 13 are for smooth stainless steel conden- 
Sees at /00, 1400, and 2800 RPM respectively. For both stain- 
Mess steel and copper smooth condensers, the following general 
observations apply: a) For the same external heat transfer 
Meer tcient, the heat transfer rate for the cylindrical smooth 
Condenser is less than the equivalent tapered condenser. 5b) As 


the external heat transfer coefficient increases, this difference 





in heat transfer rate also increases. c) The rotational speed 
Mase a greater effect on the tapered condenser heat transfer 
rate. For example, the maximum heat transfer rate of a copper 
tapered condenser will increase by a factor of 1.67 when the 
Potational speed is increased from 700 to 2800 RPM. In the 
Seandrical copper condenser, for the same change in rotational 
speed, the heat transfer rate only increases by a factor of 1.51. 

These same observations hold true for the smooth stainless 
Steel condensers. But, due to a greater thermal resistance in 
the wall of the stainless steel condenser, the heat transfer 
rates are less for all cases considered. It should be noted 
that the thermal conductivity of stainless steel is only 4% the 
thermal conductivity of copper. This accounts for the increased 
thermal resistance. 

Figures 14, 15, and 16 compare axially finned cylindrical 
With axially finned tapered copper condensers at 700, 1400 and 
2800 RPM respectively. Note that the heat transfer rates of 
the cylindrical condensers are only slightly less than those of 
the tapered condensers. This is because the heat is primarily 
transferred through the extended surface, i.e., the fin. Thus 
the film condensate in the trough has less effect on the heat 
transfer rate. In fact, the average difference in heat transfer 
rate for 700 RPM and an external heat transfer coefficient of 
Mumbtu/hr-ft*-F is 12.65%. This difference increases to 15.6% 
as the heat transfer coefficient is increased to 1000 Btu/hr- 


ft*-F. However, as the rotational speed is increased to 2800 
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RPM, the corresponding average differences in heat transfer 
Petes decrease to 12.58% and 13.35% respectively. 

Figures 17, 18, and 19 compare the heat transfer rates of 
axially finned stailness steel cylindrical condensers with 
those of axially finned stainless steel tapered condensers 
je, /00, 1400 and 2800 RPM respectively. Note the difference 
mieieat transfer rate increases as the external heat transfer 
Se@erticient increases more than in the case of the copper 
condensers. At the low heat transfer coefficient, the limit- 
Dee 


However, as the external heat transfer coefficient is increased, 


ing thermal resistance is that of the external surface IR oe 
the limiting thermal resistance becomes that of the wall due 
to the low thermal conductivity of stainless steel. 

Another observation to be noted is the fact that the rota- 
tional speed has a greater effect on the heat transfer rate of 
the cylindrical condensers than on the tapered condensers. As 
the rotational speed increases, the film thickness decreases 
due to the greater centrifugal force exerted on the film. For 
the cylindrical condenser, the maximum film thickness is much 
greater than for a tapered condenser. For example, an axially 
finned stainless steel cylindrical condenser rotating at 1400 
RPM with an external heat transfer coefficient of 1000 Btu/hr- 
ee--F has a maximum film thickness twice that of a corresponding 
eeperea COndenser. This being the case, higher rotational 
speeds will have a greater effect on the greater film thickness 


and the difference in heat transfer rates will decrease. 
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Paamese 20, 21, and 22 compare the heat transfer rates of 
copper smooth cylindrical condensers with copper axially finned 
cylindrical condensers at 700, 1400 and 2800 RPM respectively. 
imewrtigures indicate that for low external heat transfer coef- 
ficients, little is to be gained by the addition of axial fins. 
But, as the external heat transfer coefficient increases, the 
advantage becomes significant. As an example, consider Figure 
21. For h=100 Btu/hr-ft?-F, axial finning increase the heat 
transfer rate by 16%. But, for an external heat transfer coef- 
ficient of 1000 Btu/hr-ft*-F, the heat transfer rate increases 
by 194%. 

Note also, that as rotational speed increases, the advan- 
tage to be gained by axial finning decreases. This can be 
explained by the fact that for an axially finned condenser, 
the majority of the heat is transferred by the fin surface. 

As the rotational speed increases, the film thickness in the 
trough decreases. This in turn will expose slightly more fin 
surface area. On the other hand, for the smooth condenser, 
teme decrease in the film thickness will have a greater effect 
on heat transfer rate in that the thermal resistance of the 
film has decreased. Comparing the two geometries, the change 
in Overall thermal resistances for the smooth condenser will 
be greater than that for the axially finned condenser account- 
ing for the slight decrease in advantage with increasing 


rotational speed. 
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megures 25, 24, and 25 correspomd to Figures 20, 21, and 
22 but for stainless steel condensers. The results are similar 
mame COpper Situation discussed above, by the effect caused 
Meeetne increasing external heat transfer coefficient is not 
as dominant due to the high thermal resistance of the stainless 
steel wall material. 

Figure 26 compares a smooth cylindrical copper condenser 
with a smooth cylindrical stainless steel condenser at 1400 
RPM. As to be expected, the heat transfer rate of the copper 
condenser is greater than that of the stainless steel condenser 
m@ue tO the difference in the thermal conductivity of the two 
materials. The difference in the heat transfer rate is least 
for a low external heat transfer coefficient where the external 
Bi@ermal resistance 1S dominant. As the heat transfer coeffi- 
cient increases, the thermal resistance of the wall of the 
condenser becomes more important resulting in an increasing 
difference between the two condensers. 

Figure 27 provides a comparison of axially finned cylindrical 
copper and stainless steel condensers at 1400 RPM. The advantage 
of copper over stainless steel is obvious. Note that the heat 
transfer rate for the copper condenser at 500 Btu/hr-ft?*-F is 
nearly identical to the stainless steel condenser heat transfer 
fees at 1000 Btu/hr/-ft*-F indicating the advantage of copper 
Over stainless steel. 

The final analysis compared the heat transfer rate of an 


axially finned copper condenser with a triangular fin profile 
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Momtne heat transfer rate of an axially finned condenser with 
mereccangular fin profile. The rectangular fin was assumed to 
have an adiabatic tip. This comparison was accomplished for 
both cylindrical and tapered condensers. Table IV lists the 
parameters used in the analysis. These parameters are in addi- 
mom tO those listed in Table I. This comparison was only 
accomplished for the one set of operating parameters listed in 
Table IV. Note that the operating parameters chosen were the 


median values. 


TABLE IV 


List of Parameters Used in 
Rectangular/Triangular Fin Profile Comparison 


Heat Transfer Coefficient = 500 Btu/hr-ft?-F 


1400 RPM 


Rotational Speed 


Saturation Temperature 1 Z0m Ur dec recs. P 
Fin Height = 0.05906 inches 


maobe V lists the results of the analyses. 


TABLE V 


Results OL enectameular/Ilriangular 
Fin Profile Comparison 


CONDENSER GEOMETRY PLONE RORDOLE © {8tu/hr) 


tapered En Tamnole Oko S95 
eeecangile 6143.03 

Sy liner ied i triangle 5358.4 
Ge Gtamo ie 550500 
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Note that the heat transfer rate varied by only 0.4% for both 
the tapered and cylindrical condensers. The reason for this 
Meotoniiicant Varlation is provided in Table VI and VII. 

Table VI lists the convective surface temperature distribution 
eaea symmetric unit section of a tapered condenser at the middle 
increment of the condenser. Temperature location 1 is located 
Memene tip of the fin. For the rectangular profile, temperature 
meeation 1 is located at the intersection of the adiabatic 
surface (the tip of the fin), and the vertical surface of the 
fin. Temperature location 5 is located at the base of the fin. 
Locations 6-9 are located in the trough region and the remaining 
temperature locations are along the external surface of the 
Secti0On. Thus temperatures 1-5 provide the temperature: distri- 
bution along the convective surface of the fin. Note that the 
bemperature distribution is lower for the rectangular fin. 

Thus, the driving force for heat transfer, the temperature 
difference between the saturation temperature and the surface 
mameerature is greater for the rectangular fin. Thus, in spite 
Of the fact that the fin surface area for heat transfer has de- 
creased by 11%, the average surface temperature difference has 
increased by 25%. It should also be noted that the average 

heat transfer coefficient for the fin surface for the rectan- 
gular fin also decreases by 12%. However, the increase in temper- 
ature difference is the dominant change that only allows a 0.4% 


decrease in heat transfer rate. 
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TABLE VI 


Surface Temperature Distribution for Rectangular 
anamiubanemtar Axtaiiy Pinned Tapered Copper 
Condensers at Middle Increment of Condensers 


Temperature Triangular Fin Rectangular Fin 
focation Temperature (F) Temperature (F) 
i 1Lises5 Mey co's 
2 118.45 ives s 
3 118.04 Laos 
4 Lee) Teo 
5 7 705 ro 59 
6 PE Ozo 2 HHO 4:5 
i Gry G1 6). 35 
8 EPG 359 Ou 5 
9 PO .4S E16 53754 
EO IG 7, ieOn oS 
ie 116.46 ESL) ape) 
if 2 116.46 Helee 2 
ES 116.44 eb O.. 2.0 
14 116.42 Oe 57 
5 Gro. leone Z 
16 lG 3.5 1 Ory 
7 MGs 27 116.04 
18 TGA 26 LeeG0:S 
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TAS iasy It 


Surface Temperature Distribution for Rectangular 
MiGuintanodlar Axtally Pinned Cylindrical Copper 
Condensers at Middle Increment of Condensers 


Temperature lite eae Claes fem Rie Greco ued ly it Ti 
bocation Temperature (F) Temperature (F) 
1 Ieore 3.5 Tali, ore 
2 aoe 0S 117.44 
i 7 eo AO ages 5 
4 TAG? pele ore 7 
5 ila 5 1 G=25 6 
6 Ones LG 
if Gees 6 207 
8 Dele eyez dao 0M 
9 oe. 20 115.99 
10 eae 0 Ss L550 
Tea eG OZ ia 29 
EZ Pron UT ioe 9 
iS IS eS) Liat! 
14 ese 97 se 5 
ES tis .04 Saal 2 
16 IS) Sah lesieco? 
7 ieieorgcrs dees 0 7 
18 iso i500 
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liapre VIL andi¢ates a similar Situation exists in the 
mmeeidrical Condenser. Again mote the decrease in the average 
surface temperature of the fin for the rectangular profile. 
The surface area and heat transfer coefficient has decreased, 
But the increased temperature difference compensates for these 


changes, limiting the overall decrease in heat transfer rate. 
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Figure 12. Heat Transfer Rate versus Saturation Temperature 
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Figure 13. Heat Transfer Rate versus Saturation Temperature 
HeoncGmerapened andeGyrindrical Stainless 
steel Condensers at 2800 RPM 
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Figure 14. Heat Transfer Rate versus Saturation Temperature 
POGextal iy finned Cylindrical and Tapered Copper 
Condensers at 700 RPM 
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Figure 15. Heat Transfer Rate versus Saturation Temperature 
Prom xalivwrtinned Cylindrical and Tapered Copper 
Condensers at 1400 RPM 
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Figure 17. Heat Transfer Rate versus Saturation Temperature 
for Axially Finned Cylindrical and Tapered Stain- 
less Steel Condensers at 700 RPM 
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for Smooth and Axially Finned Stainless Steel 
Cylindrical Condensers at 2800 RPM 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

itiesheat transfer rate of a cylindrical condenser is less 
than an equivalent tapered condenser. This decrease in heat 
transfer rate is most Significant in a smooth condenser, where, 
femending on the external heat transfer coefficient and rota- 
tional speed, can be as great as 45%. This decrease in heat 
Sransiter rate becomes less Significant for an axially finned 
memdenser where the average decrease, fcr the range of heat 
transfer coefficients examined, was 13% for a copper condenser 
fade 1st for a stainless steel condenser. When such factors 
as the cost and difficulty in manufacturing tapered axially 
finned condensers are considered, this 13% decrease in heat 
transfer rate becomes tolerable. From a practical standpoint, 
development and analysis of cylindrical axially finned con- 
densers should be encouraged by these results. If environ- 
mental conditions permit, copper should be preferred over 
eeeinless steel due to its exceptionally high thermal con- 


ductivity and resulting higher heat transfer rate. 


B. RECOMMENDATIONS 
1) Build and experimentally test both smooth and axially 
finned cylindrical condensers to obtain experimental data 


ror comparison with results of this analysis. 






2) Develop models for rectangular and trapezoidal fin 
Meet leS With non-adiabatic tips and incorporate into code. 

3) Experimentally test axially finned cylindrical con- 
densers with rectangular and trapezoidal fin profiles to 


Sbtain data for comparison with theoretical results. 
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APPENDIX A 


Powe nOPlLLCESEENi TE ELEMENT SOLUTION 


A. SMOOTH CONDENSER 
ite analysis of the film profile in a smooth cylindrical 
condenser developed in Chapter II resulted in the following 


ordinary, nonlinear, second order, differential equation: 


d ae a 3ke(Teat - Tyw)u 
x * =e sa Ww 
es fra gs] ar (eqn A.1) 
Pwr he, 


This equation can be rearranged and expanded to yield: 


2 


* 
go ik (cane 25.2) 


dace dx 
where kK = Be eee - fu 


eae Re, 


The statement of the problem for the formulation of the 


Finite Element Method is: 


5%” d*6* fe 
ade 


. 2 
Soo: (ge) = -K (eqn A.3) 


with the following boundary conditions: 


Pymeat x = 0, ¢* = ee 


I) 


Eee at x = 0, dsg*/dx = 0 
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iicmc@onaime (length of the Condenser) is divided into ele- 
ments Of length Ax with the exception of the first and final 
elements which have length Ax/2, where Ax is the length of the 
condenser divided by the number of axial increments (NDIV). 
A system nodal point is located at each end of an element. 
Thus, system nodal point 1 is located at x=0 and the last sys- 
tem nodal point, which is equal to the number of elements plus 
1 is located at the overfall into the evaporator, x=L. All 
Meecinal nodal points are located at a position corresponding 
Menthe midpoint of the axial increment. 

Define the approximate value of the film thickness in the 


following manner: 


h 

a * 

5 (ec =") Gad. = Gd (eqn A.4) 
I 


where 6 = the approximate value of the film thickness (6*). 
G; = sene Globally basis £unctions. 
n = the number of system nodal points. 


gece the solution vector. 


On an element level, equation (A.4) becomes: 


ad. tean, Ass) 


where o> CieCelOcaibasis TUnEtions. 
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Define two degrees of freedom at each nodal point, i.e., 
both 6 and dé/dx are continuous. Thus, the local basis func- 


tions used in the finite element solution are: 

















ee Saees 
8] t g2 : Q3 (eqn A.6) 
2 2 3 
Pe eo (eqn A.7) 
ib Gee 
See 2G 
gz = iS - i. (eqn A.8) 
2 3 
a= Sa a ae (eqn A.9) 
x Q 
where & = the length of an element. 
g, and gz = magnitude Dass Lune t LOnSr. 


g, and g, = slope basis functions. 
Substituting the approximate film thickness 8. into the 


differential equation ({A.3) results in: 


4 2 ss 2 
6 d"Se 433 aoe ee (eqn A.10) 
e e e 
dx dx 





To remove the "nonlinearity" from the problem, equation (A.10) 


1s modified inthe following manner: 





+ 3n?n' ——= -K (eqn A.11) 


ag 






n is defined as the approximate value of the film thickness 
from the previous iteration. In like manner, n‘' is defined 
as the approximate value of the rate of change of film thick- 
meso with respect to x from the previous iteration. 

Forming the residual of equation (A.11) yields: 


T 


R, = n*(G' a) + 3n’n'(G) d) + K (eqn A.12) 


Invoking the Galerkin criterion for the determination of the 


Bemution vector d, l.e. 


fG;R dx ie le 5, an (eqn A.13) 
yields: 
4 i. a Te = 
faeeGG ddx + Ssn-n'f GG ddx + KsfGdx = 0 (eqn A.14) 


Each of the integrals in equation (A.14) are defined in the 


following manner: 


n 
roc ddx = 2 gg d_ dx (eqn A.15) 
x i. il 


fgg d dx = f 83 , d dx (eqn A.16) 


This integration results in the following 4x4 local elemental 


A matrix for any element: 
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and on an elemental level: 


i 
Jygg d dx 


Ly 


an 
2) 


a 


' 
On 


J- 4 


<8] g5 gz gy > 


ioe 


rb] 
4 


iE 


d 
e 


Gine 


(eqn A.17) 


ec) Awol 3) 


Gimanintegration results in the following 4x4 local elemental 


Beatrix for any element: 


- 


e 
“IU 
0 


a, 2! 


4 
es Gp 2 


de de de 


j 
+ 


ae 


I 
wR 






Mmesely, let 


n 
£Gdx = 2 I, gdx Reqgim A. 19) 
7 = 


or, on an elemental level this becomes: 


$1 
2 

Lagdx = Lo 3° dx fecame. 20) 
4 


ies integration results in the following local elemental 


eouumn F vector: 


ba[z> 


;4 
ms, 


a 


Thus, for a given element, equation (A.12) becomes: 


uy 7 3 $ = . 
Ee PY a) | ae [By d. ae iene A521) 
As mentioned above, n and n' are the approximate values of 6¢* 
and dé*/dx respectively from the previous iteration. For the 
miitial iteration, n is set equal to Oe and n‘' is set equal 

Eo 0. 


Each elemental matrix is placed in a global system [A] 


matrix with the location in the matrix based upon the local/global 
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meee point correspondence. For example, for element Nr. 2 

moe nodal points corresponding to local nodal points 1,2,3, 

aaa 4 are global nodal points 3,4,5, and 6 respectively. 
Therefore, the sum of the two 4x4 elemental matrices ([A] and 
[B]) and multiplying constants will form a single 4x4 local 
Matrix. Element (1,1) of this local matrix will be placed in 
the global A matrix location (3,3) and element (4,4) of the local 
elemental matrix will be placed in the global A matrix location 
(6,6). Ultimately, the following global system would be assembled: 


[A] on Se (eqn A.22) 


MXM 


Note that the global A matrix would be an mxm size matrix 
where m is equal to twice the number of system nodal points to 
mememic for the two degrees of freedom at each nodal point. 
Semaesoimilar fashion, d and F would be column vectors of size m. 

Once the system is assembled the boundary conditions are 
applied. This is done in the following manner: For boundary 
Sengition (a), A(1,1) is set equal to 1.0 and the remaining 
emeementS in the first row, i.e. A (1,1) 1=#=2,3,4,...m are set 
equal to 0.0. Then F(1) is set equal to Opes See is ini- 
tially determined by a relationship developed by Leppert and 
Nimmo [Refs. 8 and 9]. This establishes the value of d(l) as 
ee meeli a similar manner, boundary condition (b) is applied 


by setting A(2,2) equal to 1.0 and all other second row elements 


of the global A matrix are set equal to 0.0. Then, F(2) is set 


eal: 






Pioieto 0.0, his establishes the solution vector d{Z) as 0.0. 
Meme that d{Z) corresponds to the slope at the first nodal point. 
One additional boundary condition is required; this is the 

value of the film thickness at x=L. This boundary condition is 
necessary to completely specify the problem. The value of the 
film thickness at the overfall may take on any value depending 

on the geometry at the overfall. In the case of this analysis, 
this value was taken as 0.25 - aan This value was chosen for 
the following reason: Leppert and Nimmo [Refs. 8 and 9], in 
feoimilar analysis for laminar film condensation on a horizontal 
Surface derived an analytical solution to equation (A.1), assum- 
imeoma CONStant surface temperature. They found the film profiles 
for 6* with the overfall value less than 0.40 .- nea were essen- 
tially constant and thus any value of §* at the overfall less 
than 0.40 - ee would result-—in the same profile. In the veri- 
meeati0On Of the finite element solution, not only was this found 
to be the case, but it was also found that the heat transfer 

rate was relatively insensitive to the shape of the film profile. 
In fact, it was found that the film thickness at the overfall 
could be increased to a value as great a 0.90 - ieee and the 


resulting variation in heat transfer rate was only 10%. This 


being the case, the value of 0.25 - §* 


Was arbitrarily chosen 
max 


nOT ee 
The third boundary condition was applied by setting A(m-1,i)= 
0.0 where 1=1,2,3...m. Then, the global matrix element A(m-1,m-1) 


Was set equal to 1.0. Finally, F(m-1) was set equal to CO ee 


OZ 





Once all three boundary conditions were applied, the sys- 
tem given by equation (A.22) was solved for d. The values of 
miemapproximate film thickness, i.e., d(i), i=1,5,5,...m-1 are 
then compared to the values of film thickness from the previous 
iteration. If the relative difference is less than or equal 
BO a Specified convergence criterion (i.e., 0.0001) at all 
nodal points, convergence is met and the latest values of d 
are the solution values of &%. 

If convergence is not met, the values of d are saved for 
the next iteration where they are used to determine n and n' 
as discussed above. This iterative process is continued until 
convergence is met or until a maximum number of iterations have 


occurred. 


B. AXIALLY FINNED CONDENSER 

The finite element solution for the film profile in an 
axially finned cylindrical condenser is very similar to that 
Of a smooth cylindrical condenser. For this reason, only the 
Variations in the development will be addressed. From Chapter 
mee che differential equation for mass flow rate in an axially 
finned cylindrical condenser is given by: 


d Gea 


5° ae is (e6*°+ s*'tana)] = - Sk¢(Tsat - Tw) ue 


AD 
PO, wr Neg 
Ake ( T = jf Naz 1 
Mecosg | Sat avg" (eqn A.23) 


OP, wr he cosa 


g 


OS 






This equation can be rearranged and expanded to yield: 


d*6%* si ae alcn: acon: x4 ig 
=a (ed*"+ 6 Ue eo ae tO wand =) 





a. x 
= -K, -K,6 (eqn A.24) 
where K = Ket oat - THE 
0. wr Neg 
ahs 
4k -(T Tayo) HZ 


K, = ZEOS _ i sat avg 
0. wr hee cosa 


Supstituting equation (A.5) into equation (A.24) results in: 











eo. 06} dé. cn ' do, 
72 (ed, + 5, tana) + a (3e6 a on tana a) 
= Ky , K5o, (eqn A.25) 


Memeremove the "nonlinearity" from the problem, equation (A.25) 


is modified in the following manner" 


ds . ds. 
(ey*+y°tana) + 
dx? dx 








(je, 7 4y ¥ Panc) 


= - kK, - Ky (eqn A.26) 
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Here, y and y' are defined by the following relationships: 


m= Or z R*(6, -65_4) (eqn A.27) 
yt = 3, + R*(S, -3,_,) (eqn A.28) 
where 6. = approximate value of the film thickness for 
the present iteration. 
eo approximate value of film thickness from the 
DRev lous 1te ration. 
ce = approximate value of the derivative of the film 
thickness for the present iteration 
Oa ip Gowiidtemvalile Gi the derivative Of the £1 1m 
thickness from the previous iteration. 
hee= relaxation factor. 


These two variables, y and y', are in actuality, adjusted approx- 
mation of film thickness and derivative of film thickness re- 
Spectively. This adjustment is required in order to converge to 
a solution. 

The finite element solution is now identical to that of the 
smooth condenser, that is, the residual is formed, the Galerkin 
criterion is invoked, and identical 4x4 local matrices are 


Memived. Finally, the equivalent of equation (A.21) is formed. 
4 5 Sat iat ° 
Meee Y, tance) [Aj], + (Sey,°¥, + Y, Y, tan) [Bl] °4, 


a ( - Ki 6 - Ky a) iE (eqn A.29) 
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Notice that this equation has two forcing terms. The additional 
term (K5Y) resulted from the nonlinearization of the problem. 
wust ai in the smooth condenser film profile solution, the 
elobal system given by equation (A.22) is formed, the boundary 
conditions applied and the system solved for a solution vector 
Mueitee tterative process is continued until convergence is met. 
With each iteration, y and y' are updated and used for the 
next iteration. When convergence is met, the latest values of 


Mee, 2=1,3,5,...m-1 are the solution values of 6*. 
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APPENDIX B 


USER'S MANUAL 


This appendix describes the data cards required to use 
the computer code. 

imewdata 1s divided into “blocks" for convenience. Each 
page of this user's guide is a separate block. For each block, 
a general description, the required format, and appropriate 
comments are provided. 

It is imperative that input data be consistent or errors 
will result. For example, if a smooth geometry is being analyzed, 
the finite element parameters must also result in a smooth model. 
In addition, all data fields must be filled with an input value, 
Perene if that value is not needed for the analysis. For example, 
in a smooth analysis, no fin half angle is required for the cal- 
Gulations; however, a value of the correct format must be pro- 
vided in the fin half angle field or an INPUT/OUTPUT error will 


result. 
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DATA BLOCK A 


MescRIPTION: FINITE ELEMENT PARAMETERS 


FORMAT : 


—_— —- @=& =e =e =| = 


_-_ = = = = = — 





8I5 


-_— se =e we we wees ws wSse Sse Sse ws BP |= ws Sse we wee#® fs | — BSF SF BS = Sse = = SB BF BF BS BS fF 2 £2 =—§ = |" ese ee se ee Pe we ese fF fF 2s |S - -— 


NCMREC NCMTRI NRWFIN NRWTRF NCMTRFE NCOL NPRNT 


nae oem el rll rl il rl rll rl eclectic i lcrrhlClC rerclUlClcrecmlUClcrechlClcrermlUlClcr reemlCl llc rel cell eel eet ll rece ee ll ell ect eel ell ret eee rel etl ell rschlcr rer lr Pr lc rl ell el Eh 


= = == eee eee ee eee ee ee ee 


CONTENTS 


NDIV---Number of axial increments. Must be less than 
O. equal to 50. 


Cem columns Of tinzrte €lements in the rec- 
Panganiban SE€Ction of fin. May be “equal to 0 if 
fin 1s triangular only. 


Neila Nninbes. o2 columns of finite elements in tri- 
ANCGeareseehr Ono, tim. May be equal to 0 if 
OMiveereetangular ot Mn. 


PemeiN-Numbem of rows Of finite elements in the fin 
section of model. May be equal to 0 if a smooth 
Gonde tise 1. 
Noe we lieerlangulear OG trapezoidal fin, NRWFIN 
must equal NCMTRI. 


NRWTRF-Number of rows of finite elements in wall section 
of the model. 


NCMTRF-Number of columns of finite elements in the 
trough section of finned model. Set equal to 
O 1f a smooth condenser. 


NCOL---Total number of columns of finite elements. 
Must be equal to NCMREC+NCMTRI+NCMTRF for a 
finned model. 


NPRNT--Print control number. 
PwcdudlmcoulmcorGespondence table and major 
elements Of finite element model will be printed 
eLig 
If equal to 0, output will be suppressed. 
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DATA BLOCK B 
DESCRIPTION: FLUID, FIN MATERIAL SELECTION PARAMETERS 


FORMAT: 215 


EE LD CONTENTS 


1 MeeUID-It equal to Ogeworking fluid is water. 
If equal to 1, working fluid is freon. 

2 IFIN---If equal to 0, condenser wall material is 
COpper. 


If equal to 1, condenser wall material is 
stainless steel. 


OS 






DATA BLOCK C 
DESCRIPTION: CONDENSER GEOMETRY 


FORMAT: 6G10.5 


_—- =e = ee we ww we wT BF BSF WM" Bs BF BS BS K|— Fe Fe S| S| Ss SF Se Bese |S — SF ww KX Pe SF ww Se S| Se Sse ew FF |B BF FB FF T_T Ft Fe Te eT cell cel eS OO el ee ll rel ee OO ele hl 


CLI Rewoel THICKI Be NL CANGL FNWTHI 


en i i i i ee ed 


Ere LD CONTENTS 
1 moe --condenser length, (inches). 


2 RBASEI-Inside radius to wall of condenser at condenser 
end (inches). 


3 THICKI-Wall thickness (inches). 


4 BFINI--Fin height (inches). Must be set equal to 
0.0 if smooth condenser. 


3 CANGL--Condenser half angle for tapered condenser 
(degrees). Must be set equal to 0.0 for 
cylindrical condenser. 


6 PNWIHI-Width of rectangular portion of trapezoidal 


Hepecehanieuhateatin (imches). lf triangular 
tiny set equal to 0.0. 
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DATA BLOCK D 
DESCRIPTION: INTERNAL FIN GEOMETRY 


BeRMAT: 2G610.5 


= = = wee fF =e @e =e we ee se Ss |e |S =e ee |e Pe | | Be S| |e PF Ss Ss |S Ss Pw Pe et ll reel ree ll eet lll eel etl ee ell etl ee ll ee lll ee ll ee ll ell ee ee ee eel ee ese Pell Or ll rll ee ll ee ll esc ll ell SS eS el 


FANGL BoE © 
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FLELD CONTENTS 


i FANGL--Fin half angle (degrees). Set equal to 0.0 
for smooth condenser or rectangular fin. 


2 ELOEO-=Retionot seroqugh width to fin base width. 
Determines spacing between fins. 






DATA BLOCK iE 
DMipoeRIPTION: CONVERGENCE CRITERIA 
meeeAT: 2610.5 


—=— = se we we |B |B BO VB! Ss | SF SF SB ks —Bs | Bs | Hs ws Ss SF SF |B Ss Kw | | we Sse BB Sse B |B se BF BS BS BS BS Sf = SF fF fS Sf |] SF ff] |] = SF SF Ss | SF fs =F fF 2-2 —- — 


_se we =e ws = SB ss BS | | & | SF SF SF S| we BS we SF | SB |] |S BS VBS] BS Bs BSB BSB | Fre Fs SF Ses Fe BF SF SF SF w_eZseee swe SF SB SF Se Se = SF SF SF SF SF SS | Ss — 


— — — = ee eee iii eee sl eee le le el el el ell elle el ell el el ll el 


Bee 1D CONTENTS 


1 ew --“hemm@erature convergence criterion. Used to 
determine solution of two dimensional steady 
Suamemicat COnNductIOn probiem. 


2 CRITDL-Mass flow convergence criterion. Used only 
in cylindrical condenser analysis for mass 
Plow convergence fest - 






DATA BLOCK J 
PESCRIPTION: OPERATING PARAMETERS 


meRMAT: 5610.5 


= =e wee we = Sse Sse es =F BP se sFwexexee= ww8 wBeewrew=_eew=®=ew== Bereeewe=wse 2a =e wswreese 2 w@=@ewr= =e esweewreeewe=e wewewewese @B |=& wo 2B we we ]=# ese =F @2ee 2 2S Be est ese 2s |= = = 


RPM HINF TINTL TSAT EN 


ee ee a el 


eee LD CONTENTS 
1 RPM- ---Rotational speed {revolutions per minute). 


2 HINF---External heat transfer coefficient 
(Btu/hr-ft -deg F). 


e Pel ite letempematuame estimate (degrees F) . 
4 TSAT---Saturation temperature (degrees F). 
5 TINF---Ambient temperature (degrees F). 






DATA BLOCK G 
DESCRIPTION: OUTPUT PRINT CONTROL 


FORMAT: 415 


—=—_se =e we wes =e Be Be Bese ws ws ws we wes wees se swe wese Fs =sawese |= BSB wae Fe =F se Bs ws se sews se se |= ses = F= = = se | ee ese = =F |e fF SF S| SSeS lel 


IUNITS NFLAGL NFLAG2 NFLAG3 


—-—=_= = ese fF =e se ese = ws BP wee Sse Ss BS |= ws fs BS |] Ba we ee BS Ss BFS 2S 2S 2S 2S B= 2S Fs 2S BS BS =F fF BFS TFT se se — Be we H—| ese |= =e wes ese ses ws ww se | |] = 


Eee LD CONTENTS 


iL tis -Oumput units control number: 
ceo Moe nOomcalculated results will be pro- 
Wwidcamim English units . 
Coo Ma tnput parameters will be repeated 
in SI units and calculated results will be pro- 
Vaided In ol Unlts. 
EUs = e242, nnutewanameters will be repeated 
f(Pnices wand OUlpUE results will be provided 
meoeGtmeeno lish and SI units. 


2 NETAGI- The Stirst, axial imerement at which the parameters 
PcotCamundeneGeilarksewilll be provided aS output. 
3 NFLAGZ-The final increment at which the parameters 
listed under remarks will be provided. 
4 NFLAGS-The step change in increments between NFLAG1. 
REMARKS 


1) No matter what value of IUNITS is used, input parameters 
Waetl always first appear in English units. 

2) For increments indicated by NFLAG values, the following 
parameters will appear: a) x-coordinate, y-coordinate and 
temperature at each nodal point, b) length of element, heat 
mmansrer coefficient and heat rate per unit length for each 
convective boundary element. 

3) As a minimum, the values of 1, 1, 1 must be provided as 
input for NFLAG1, NFLAG2, and NFLAG3 respectively. 
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watA BLOCK H 
MeockIPTION: TAPERED CONDENSER SOLUTION METHOD 


FORMAT: 115 


PLeELD CONTENTS 


1 NSOLVE-Tapered solution control number. 

For tapered, axially finned condenser, set 

NS ORME ly 

For tapered smooth condenser, NSOLVE must be 

set to one of the following three values: 

Scum oeivb =. 2 1f solution of film thickness is 
to be based on Ballback's [Ref. 1] 
equation. 

SeumoULVEe- Sif sselutron, of film thickness is 
to be based on Daniels and Al-Jumaily [Ref. 13] 
equation, neglecting drag terms. 

set NSOLVE = 4 if solution of film thickness is 
to be based on Daniel's and Al-Jumaily [Ref. 
13] equation, with drag effects included. 


REMARKS 


1) NSOLVE is only used in tapered condenser analysis. 





DATA BLOCK it 
DESCRIPTION: CYLINDRICAL CONDENSER ANALYSIS PARAMETERS 


meRMAT: 315,2G15.10 


_-_ = = ws w= =—§ 2 = 2s =F 2B BF BF Bs BF BF BF BFS |= = ws BF BS B| BFS B8| BF BS BS |= Bs |—§ fs BSF B28 fs |= BF |= 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 82 2+-° —- — 


MONCE ETERMX ITRPRT RELAX DELMAX 


nn omnia ae ae Ss ee ae 2s we = 2 2 2 2 =F = |= =| 2 2s |= = |= |= |= 2 |S | |S |S |S |S |S |e |]! | |S = B= |S | |= BF Bf |] = B|s |—§ S&F = Vf fF fs |] = |—§ = = = =—- —_— — 


Ee LD CONTENTS 


i NWONeGr--simele mleration parameter. 
If NONCE = 1, only one iteration will be 
permitted. 


Pee NOMGEs= 0) 1terations Wild continue until 
convergence or maximum number of iterations 
is reached. 


2 ITERMX-Maximum number of iterations permitted in analysis. 
5 MVE tweeehatlOn print Contnel parameter . 
if ITRPRT = 1, mass flow convergence test results 
will be provided for each iteration. 
If ITRPRT = 0, mass flow convergence test results 
Will only be provided on final iteration. 
4 RELAX--Relaxation variable used in finite element solution 
Wmeevlindrical finmed film profile- 
5 DELMAX-Initial estimate of maximum film thickness used in 
Soler on Gf cylindryveal finned film profile. 
REMARKS 
1) Above parameters are only used in cylindrical condenser 
analysis. 
2) Recommended value of RELAX is 0.80. It is sometimes necessary 


Bomdajust this value plus or minus 0.05 to reach film profile 
convergence at small film thickness values. 

3) Input value of DELMAX is only used in cylindrical finned 
analysis. For the cylindrical smooth condenser, DELMAX is 
internally generated. 
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